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ABSTRACT 


This  study  considers  the  self-consistent  modification  of  the  trapping 
dynamics  of  a  low  density  cold  electron  beam  du  to  the  external  appli 
cation  of  a  DC  electric  field.  For  DC  fields  smaller  than  the  peak 
amplitude  of  the  saturated  beam-plasma  Instability,  the  beam  energy 
remains  clamped  while  the  wave  amplitude  grows  secularly.  Large 
travelling  potential  wells  appear  and  craete  strongly  focused  charge 
clumps.  By  considering  the  role  of  wave  dissipation  an  exact  dynamic 
B6L  (Bernsteln-Green-Kruskal )  equilibrium  Is  found  analytically.  It 
consists  of  a  singular  charge  clump  which  propagates  through  the  me¬ 
dium  at  constant  velocity  even  though  a  DC  field  Is  present.  The 
numerical  study  of  Ihe basic  equations  shows  that  the  system  evolves 
asymptotically  to  this  singular  state.  A  variety  of  experimentally 
relevant  phenomena  associated  with  the  trapping  dynamics  Is  investi¬ 
gated,  Including  the  stability  of  the  dynamic  nonlinear  equilibrium 
to  sideband  growth. 
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This  study  considers  the  self-consistent  modification  of  the  trapping 
dynamics  of  a  low  density  cold  electron  beam  due  to  the  external  application 
of  a  DC  electric  field.  For  DC  fields  smaller  than  the  peak  amplitude  of  the 
saturated  beam-plasma  instability,  the  beam  energy  remains  clamped  while  the 
wave  amplitude  grows  secularly.  Large  travelling  potential  wells  appear  and 
create  strongly  focused  charge  clumps.  By  considering  the  role  of  wave  dis¬ 
sipation  an  exact  dynamic  BGK  (Bemstein-Greene-Kruskal)  equilibrium  is 
found  analytically.  It  consists  of  a  singular  charge  clump  which  propagates 
through  the  medium  at  constant  velocity  even  though  a  DC  field  is  present. 

The  numerical  study  of  the  basic  equations  shows  that  the  system  evolves 
asymptotically  to  this  singular  state.  A  variety  of  experimentally  relevant 
phenomena  associated  with  the  trapping  dynamics  is  investigated,  including 
the  stability  of  the  dynamic  nonlinear  equilibrium  to  sideband  growth. 
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I .  Introduction 
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The  trapping  of  charged  particles  within  the  travelling  potential  wells  of 

collective  plasma  oscillations  is  one  of  the  more  significant  strong  nonlineari- 

1  2 

ties  encountered  in  plasmas.  Particle  trapping  is  known,  both  theoretically  ’ 

t 

and  experimentally  ,  to  alter  the  linear  collisionless  damping  (Landau  damping) 
of  collective  modes  supported  by  plasmas,  and  it  is  one  of  the  dominant  mechan¬ 
isms  that  limits  the  growth  of  velocity-space  instabilities.  The  strong  effects 
produced  by  particle  trapping  are  related  to  the  formation  of  phase-space  granu¬ 
lations  out  of  initially  smooth  contours.  The  phase-space  granulation  is  par¬ 
tially  reversible,  and  under  carefully  selected  conditions  can  lead  to  inter- 
4  S  S  6 

esting  amplitude  *  and  phase  '  oscillations  which  persist  for  relative  long 
times  of  experimental  interest.  The  phase-space  granulation  can  disappear 
smoothly  by  the  process  of  phase  mixing  while  the  charged  particles  remain 
confined  within  the  potential  wells,  as  is  characteristic  of  the  time  asymptotic 
formation  of  steady  Bemstein-Greene-Kruskal7  (BKG)  modes.  The  granulation  can 
also  disappear  more  abruptly  if  the  trapped  particles  escape  from  the  potential 
wells,  as  may  occur  due  to  the  spontaneous  growth  of  sideband  waves,  sudden 

g 

phase  slippage  due  to  finite  dissipation  ,  or  the  launching  of  other  large 
o 

amplitude  waves. 

The  present  study  is  concerned  with  the  dynamics  of  trapped  particles  and 

the  associated  phase-space  granulation  in  the  presence  of  an  externally  applied 

DC  electric  field.  The  motivation  for  investigating  this  dynamical  system  in 

detail  can  be  attributed  to  the  finding  in  a  previous  theoretical  study10  that 

the  momentum  of  a  runaway  beam  can  become  clamped  due  to  the  excitation  of  modes 

travelling  in  synchronism  with  the  beam.  That  earlier  study  was  in  turn  stimu- 

11  12 

lated  by  the  experimental  observation  '  of  runaway  energy  clamping  in  the 


Microtor  Tokamak  at  the  University  of  California,  Los  Angeles. 
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As  is  expected,  the  beam  clasping  effect  sets  in  When  the  amplitude  of  the 
resonant  wave  becomes  larger  than  the  external  DC  electric  field.  The  earlier 
calculation  of  this  effect  has  been  accoaplished  by  aeans  of  a  spatially  averaged 
formalism  based  on  the  exact  conservation  laws  for  energy  and  momentum,  but  with 
the  restrictive  iaplicit  assumption  that  the  phase-space  granulation  is  instantan¬ 
eously  saeared  out,  i.e. ,  the  runaway  bean  acquires  an  irreversible  thermal  spread. 
Of  course,  in  reality  the  trapped  particle  dynaaics  is  partially  reversible. 
Therefore,  it  is  of  interest  to  ellucidate  the  underlying  single  particle  dynaaics 
accompanying  the  damping  process,  to  validate  the  existence  of  damping  indepen¬ 
dently  of  the  spatially  averaged  formalism,  and  to  establish  the  limitations  of 
■  .the  latter.  In  addition,  there  are  various  associated  effects  which  are  of 

interest  foT  both  basic  and  applied  reasons.  A  preliminary  study  of  some  of  these 

13  14 

effects  has  been  previously  discussed  *  by  this  author. 

To  isolate  the  principal  effects  produced  by  an  external  DC  electric  field 
Eq,  we  consider  the  cleanest  environment  for  trapped  particle  dynamics,  i.e., 
the  low  density  cold  beam  problem  in  which  the  background  component  supporting 
the  waves  remains  strictly  linear.  In  practice  this  situation  can  be  realized 
by  having  a  cold  background  plasma  with  heavy  ions  or  simply  replacing  the  cold 
plasma  by  a  slow  wave  structure,  as  is  commonly  done  in  travelling  wave  tube  (TUT) 

A 

devices.  The  cold  plasma  assumption  can  be  met  by  supplying  sufficient  neutral 
gas  in  order  to  compensate  for  the  heating  produced  by  the  waves  and  the  intrinsic 
ohmic  heating  associated  with  Eq. 

It  is  well  known®* 15’1^  that  in  the  absence  of  the  DC  electric  field  the 
weak  cold  beam  amplifies  a  narrow  spectrum  of  noise  wates  supported  by  the 
background  structure.  After  several  e-foldings  the  noise  evolves  into  a  nearly 
monochromatic  signal  (single  mode)  whose  growth  is  Stepped  by  the  abrupt  trapping 


of  the  bean  electrons  within  its  potential  troughs.  The  application  of  a  DC 
electric  field  to  this  system  is  found  to  produce  a  secular  growth  in  the  wave 
amplitude  without  a  noticeable  increase  in  the  beam  momentum,  i.e. ,  the  external 
momentum  push  is  transferred  to  the  wave  while  the  beam  remains  clamped.  This 
result  is  obtained  in  the  present  study  by  following  the  self-consistent  orbits 
of  the  trapped  particles  in  a  digital  computer,  thus  confirming  the  existence  of 
beam  clamping  predicted  by  the  earlier  spatially  averaged  analysis.10  However, 
additional  features  associated  with  the  trapping  dynamics  are  uncovered.  One 
of  them  is  the  expected  appearance  of  trapping  oscillations  of  increasing  fre¬ 
quency  which  modulate  the  secular  growth  in  the  wave  amplitude.  Another  feature 
is  the  generation  of  secondary  low  density  runaway  beams  which  arise  from  the 
spillage  of  the  primary  beam  particles  out  of  the  potential  wells. 

In  the  absence  of  wave  damping  the  secular  increase  in  the  wave  amplitude 
produced  by  the  DC  electric  field  is  unbounded.  Large  potential  wells  are  gen¬ 
erated  and  the  trapped  particles  evolve  into  phase-space  clumps  of  ever  decreas¬ 
ing  dimensions.  When  a  wave  dissipation  channel  exists,  one  obtains  analytically 
an  exact  dynamic  BGK  equilibrium  consisting  of  a  singular  charge  clump  which 
is  being  pushed  by  Eq,  but  whose  momentum  remains  constant.  This  unusual  and 
highly  nonlinear  state  is  found  to  be  approached  asymptotically  by  the  time 
dependent  computer  solution  of  the  problem. 

The  application  of  the  DC  electric  field  permits  the  detailed  manipulation 
of  the  position  of  the  trapped  particles  relative  to  the  bottom  of  the  potential 
wells  of  the  wave.  By  changing  the  direction  of  the  DC  field,  one  can  either 
make  the  wave  grow  or  damp.  By  means  of  this  phase-space  reconstruction,  it  is 
possible  to  reverse  the  drastic  disappearance  of  the  trapped  particle  oscilla¬ 
tions  associated  with  finite  wave  dissipation. 
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The  enhanced  growth  in  the  wave  amplitude  produced  by  EQ  can  be  used  to 

iaprove  the  efficiency  of  a  TXT,  as  was  investigated  nearly  twenty  years  ago 
18 

by  Hess.  A  similar  application  may  be  found  for  the  collective  free  elect- 

19 

ron  lasers  under  present  consideration  ,  where  particle  trapping  plays  an 
important  role.  The  results  presented  here  illuminate  some  of  the  fundaarantal 
processes  which  underly  these  practical  developments,  and  define  specific  exper¬ 
iments  amenable  to  investigation  in  the  laboratory. 

In  Section  II  of  the  paper  some  of  the  basic  properties  of  trapping  behav¬ 
ior  in  the  presence  of  a  DC  electric  field  are  examined  analytically.  The 
weak  cold  beam  formalism  is  examined  in  Section  III.  The  beam  clashing  effect 
is  discussed  in  Section  IV,  and  the  dynamic  BGK  equilibrium  is  explained  in 
Section  V.  The  stability  of  the  singular  clusp  equilibrium  to  sideband  growth 
is  investigated  in  Section  VI.  Concluding  remarks  are  presented  in  Section  VII. 


II.  Modified  Trapping  Behavior 

Before  considering  the  self-consistent  evolution  of  the  wave  amplitude 
and  phase  together  with  the  trapped  particle  dynamics  in  the  presence  of  an 
external  DC  electric  field,  it  is  useful  to  develop  a  simple  description  of 
the  lowest  order  effects  produced  by  the  DC  field  on  particles  which  are  trap¬ 
ped  in  an  electrostatic  wave  of  fixed  amplitude  and  phase  velocity.  In  this 
case,  the  equation  of  motion  for  electrons  of  charge  -e,  and  mass  m,  in  the 
wave  frame  moving  with  speed  w/k  relative  to  the  laboratory  is 

-m  il  x  =  -  +  e  E0  (1) 

dX 

where  A  is  the  wave  amplitude,  k  the  wavenumber,  and  u>  the  frequency. 

It  is  evident  from  Eq  (1)  that  the  particles  move  in  a  distorted  potential 
well  whose  potential  energy  function  U  is 

U  =  -UA/k)[(lo<E,M>  (2) 

The  spatial  behavior  of  the  corresponding  potential  well  is  exhibited  in  Fig.  1 
for  the  values  EQ/A  ■  0.2,  0.5,  0.9,  over  the  interval- irdcx^ir.  It  is  seen 
that  particle  trapping  ceases  to  exist  when  E^A,  since  only  a  single  turning 
point  is  found  for  the  potential.  This  trapping  threshold  is  a  useful  comparison 
standard  for  the  studies  of  Sec  III,  where  it  is  shown  that  this  threshold  can 
be  exceeded  by  a  factor  of  2  when  the  self-consistent  amplitude  and  phase  changes 
are  considered. 

For  Eqj<0,  the  bottom  of  the  potential  well  is  shifted  from  x  =  0  to  a  location 
xg  given  by 


Xj  =  k*1  AHM'1  (£o/^ 


(3) 


and  the  corresponding  bounce  time  t^  is 

=  t4/w° 

1/2 

instead  of  the  usual  elliptic  integral  expression.  In  Eq  (4),  (eAk/m)  is 

the  bounce  frequency  for  E  =0,  U  -U(t»0)/(eA/k)  is  the  scaled  energy  of  the 

particle,  E-kx,  and  C2  refer  to  the  turning  points,  which  are  not  symmetric. 

The  lowest  order  correction  to  the  bounce  frequency  of  particles  bouncing 

near  the  bottom  of  the  distorted  well  can  be  extracted  through  the  Bogoliuvov 
20 

procedure  applied  to  the  equation 

<L  \  *./a  (s) 

Ax' 

where  5  *kx  ,  ?-y+C  ,  and  t-w.  t.  For  small  displacements  (y|«l,  Eq  (1)  can  be 
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approximated  by 

(6-/^ (.*». V>/fc  (6) 

Taking  the  solution  to  be  of  the  form 


and  eliminating  the  secular  behavior  in  lowest  other  leads  to 


A. 

ax 


*  rlT 

jae 


ccrxe  6 


Jl  Ol  =  O 

ax 

and 

i»-.  <.CA^' -  (c^V* 
ax 

a_e  =  C^vTH'-  <£/|fc) 

ax 

which  results  in  the  modified  bounce  frequency 


.« 


C3fe-  (eteA/W)K[l-^. Mfj\' - («». A? 1 


(8) 


(9) 


(10) 


for  a  particle  of  initial  displacement  xq.  It  is  seen  explicitly  from  Eq  (10) 
that  the  application  of  the  DC  electric  field  weakens  the  restoring  force  pro¬ 
vided  by  the  wave  potential,  and  furthermore  that  bouncing  ceases  to  exist  as 
Eq  approaches  A. 

The  distortion  of  the  well  produced  by  EQ/0  gives  rise  to  the  detrapping 

of  some  particles  which  for  Eq=0  would  remain  confined  within  the  wave  troughs. 

The  time  evolution  of  the  orbits  of  these  detrapped  particles  can  be  obtained 

by  solving  Eq  (1)  in  a  perturbative  manner,  i.e.,  taking  the  lowest  order  orbit 

£  to  be 
o 

+  VXxVa. 

°  (ii) 

X  =  £o/X 


The  scaled  equation  for  the  first  order  velocity  becomes 
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£ 


JUM. V  ^(o^X  v xx*/*"J 


(12) 


whose  exact  solution  can  be  expressed  in  the  terns  of  the  Si  (sine)  end  Ci 
(cosine)  integrals,  i.e. , 

= -  Wl  ( jCiX*wl-  ***». 

ce^©0 

where, 

o0=  W-  ('iwlT/ix 

?Cc}=  +  V^/uW'* 

(14) 


The  asynptotic  behavior  sets  in 
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for  z>4 


and  the  leading  expansions  are 


Sc  %  THz  - 


Cl» 


(15) 


thus  indicating  that  the  detrapped  particles  becone  runaways  which  exhibit 
decreasing  oscillations  with  a  frequency  (E0/2A) ■  (ekEQ/2n)1^2.  The 
net  effect  of  the  wave  potential  on  the  runaway  electrons  created  by  the  de¬ 
trapping  process  is  to  modify  the  free-fall  velocity  by  an  amount  |v  given  by 


tv  =.  (<VhHMrt.fi  (16) 

over  a  characteristic  time  At^lO/ug. 

The  displacement  of  the  bottom  of  the  effective  potential  given  by  Eq  (3) 
plays  a  major  role  in  the  time  evolution  of  the  wave  amplitude.  Essentially, 
the  forward  push  provided  by  EQ  overcomes  the  particle  recoil  associated  with 
momentum  transfer  to  the  wave.  Hence,  a  particle  remains  in  synchronism  with 
the  wave  and  causes  a  secular  increase  in  its  amplitude,  as  is  discussed  in 
Sec  III.  Therefore,  an  additional  consequence  of  E^O  is  that  the  depth  of 
the  trapping  well  increases  in  time  if  Eq  is  less  than  the  threshold  value  for 
the  destruction  of  trapping  (i.e.,  E^A  to  lowest  order).  An  exactly  soluble 
simplified  model  of  particle  trapping  in  this  environment  is  described  by  the 
equation 

X  *  eE0 

(17) 

where  a  measures  the  rate  of  increase  of  the  wave  amplitude.  Defining  the 
auxiliary  variables 

«  =  (*/«kf  (Me.)Uw) 


transforms  Bq  (17)  into 
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Jl.  *  *  >t*  -  \ 


(19) 


which  Is  tho  Inhomogeneous  Airy  equation  whose  solution  for  u  particle  Initially 
trapped  at  the  bottom  of  the  potential  well  (i.e.,  z  ■  dz/dri  ■  o  at  t  •  o)  is 
given  by 


Where  Ai,  Bi,  and  Gi  are  the  standard  Airy  functions.  The  constants  a  and  b 

are 

OL  — 

b  *  (  Gi.  At  -  Gi.  Ai )  /  Ai. 

(21) 


Where  the  priae  denotes  a  derivative  with  respect  to  the  argument  and  all 
functions  are  evaluated  at  n(0)  ■  (w^/o)2^3. 

In  the  asymptotic  limit  ott»l  Eq  (20)  implies  that  the  particle  position 
is  given  by 

X  »  (e./^F  yf[! 0+stf*  (22) 

In  Bq  (22)  V  is  a  phase  angle  determined  by  the  ratio  of  a  to  b,  and  F  is  the 
corresponding  amplitude,  which  is  strictly  a  function  of  (w^/a)2^3. 
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It  is  seen  from  Eq  (22)  that  the  trapped  particles  execute  rapid  oscillations 
whose  amplitude  is  proportional  to  the  displacement  of  the  bottom  of  the  well  pro¬ 
duced  by  Eq  (i.e.,  Eq  (3)).  These  oscillations  decrease  in  time  due  to  the  com¬ 
pression  produced  by  the  ever  increasing  potential  well.  It  is  this  compression 
which  is  responsible  for  the  formation  of  the  singular  charge  clumps  associated 
with  the  dynamic  BGK  equilibrium  described  in  Sec  V.  It  should  be  noted  that 
although  Eq  (22)  applies  exactly  to  a  linearly  increasing  wave  amplitude,  similar 
results  are  obtained  for  other  time  dependences,  and  can  be  calculated  asymptoti¬ 
cally  by  means  of  the  well  known  WKB  analysis. 


ft 
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III*  Weak  Cold  Boa a  Formalism 

Consider  a  fast  electron. beaa  of  initial  speed  vQ  whose  density  nfe  and 
therml  spread  are  sufficiently  snail  to  neet  the  weak  cold  beaa  criteria.16 
The  beaa  is  assuaed  to  propagate  through  a  aediun  (a  cold  plasma  of  density  nQ 
or  an  equivalent  slow  wave  structure)  whose  behavior  reaains  strictly  linear,  so 
that  this  contribution  to  the  problem  can  be  described  analytically  by  aeans 
of  a  slowly  varying  dielectric  e[k,t»(t)}.  Since  in  the  present  study  one  is 
interested  in  the  role  played  by  the  dynaaics  of  the  trapped  particles,  the 
beaa  is  discredited  into  N  charge  sheets  whose  individual  orbits  x ^ (t)  are 
followed  in  a  digital  coaputer  self-consistently  with  the  evolution  of  the  waves. 
This  procedure  is  to  be  contrasted  with  the  analytic  description  of  the  beaa 
behavior  used  in  the  spatially  averaged  formalism. 

Defining  the  scaled  position  of  particle  j  in  a  frame  moving  with  the  initial 
velocity  of  the  beaa  by  5 ^ (t)«ko(x^ (t)-vot],  where  kQ  is  the  wavenumber  of  the 
fastest  growing  linearly  unstable  wave,  yields  the  equation  of  action 

£  =  £  -  \t,  +  c  c] 

(23) 

where  the  scaled  quantities  are 

_a  =  [ 

X*=  («•»- 

A*s  E*  >  €*(«fc./<nuf)60 

(24) 
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The  physical  meaning  of  the  quantities  in  Eq  (24)  is  the  following.  is  the 
complex  amplitude  of  a  wave  having  wavenumber  k,  which  in  the  absence  of  the  beam 
has  a  laboratory  frequency  w£*wp,  where  wp  is  the  electron  plasma  frequency 

of  the  background  medium.  is  the  wave  amplitude  scaled  to  the  saturation  level 

Ep-mQ  /eko  due  to  beam  trapping  by  the  fastest  growing  wave  kQ,  when  Eq-0.  A  is 
essentially  the  growth  rate  of  the  linear  beam-plasma  instability,  and  t  is  the 
appropriately  normalized  time  variable.  5  ^  is  the  Doppler  shifted  frequency  of 
mode  kQ  in  the  frame  moving  with  speed  vq.  refers  to  the  frequency  of  mode 
kQ  in  the  absence  of  the  beam.  Finally,  is  the  real  part  of  the  dielectric, 
i.e.,  e»er*ie^,  and  (der/du)o  is  the  derivative  evaluated  at  k«k0,w«ii»®.  By 
analogy,  (ac^au)^  refers  to  the  same  quantity  evaluated  at  k  and  w?. 

The  self-consistent  time  evolution  of  A^  is  determined  from  Poisson's 
equation,  which  in  the  slow  time  scale  approximation  (0«up)  takes  the  form 

vV\ 

(25) 

in  which 

/A.V(>€T|W>K 

(26) 

Physically,  represents  the  weak  dissipation  acting  on  mode  k.  Its  origin 
is  left  unspecified  in  the  present  work,  but  can  be  conceptually  associated 
with  background  plasma  heating,  wall  heating,  radiation  leakage,  etc. 
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The  principal  dynamical  effects  associated  with  the  presence  of  the  external 
DC  electric  field  can  be  understood  fron  the  exact  conservation  relations  derivable 
from  Eqs  (23)  and  (25),  namely 


jtUT  «  i  * 


(27) 


t  ♦  t  K  lv2 

<&XM  * 

=  e(|  tv**  (*«*>] 


(28) 


where.  Re  (fa^)  refers  to  the  real  part  of  the  complex  frequency  shift  defined 

by 


(29) 


Equation  (27)  is  the  conservation  law  for  the  total  momentum  of  the  system, 
which  is  composed  of  the  beam  momentum 


P  s  s  i 
^  X  ti 


#• 


and  the  wave  swan  tun 


(30) 
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(31) 


The  total  momentum  can  be  altered  (increased  or  decreased)  by  a  nonzero  £,  and 
is  destroyed  whenever  there  exists  a  wave  dissipation  channel  (i.e. ,  r^O). 

Equation  (28)  represents  the  conservation  of  energy  principle  in  the 
frame  moving  with  speed  vQ  relative  to  the  lab.  It  states  that  the  total  energy, 
composed  of  the  mechanical  beam  energy 


(32) 


and  the  wave  energy 


(33) 


is  altered  by  the  work  done  on  the  beam  particles  byf* ,  and  is  dissipated  when 

v°- 

For  rk-0,  the  conservation  laws  take  the  simple  form 


ax 


ivO* 8 


(34) 


A.  =  6  ft 

dl  (55) 

When  €.■(),  Eq  (34)  is  satisfied  by  balancing  the  growth  of  the  nearly  mono¬ 
chromatic  wave  with  a  corresponding  decrease  in  the  beam  momentum,  i.e.,  the 
beam  recoils.  As  the  beam  recoils,  its  kinetic  energy  increases  in  the  frame 
moving  with  speed  vQ  relative  to  the  lab,  thua  Eq  (35)  is  satisfied  by  the  appear- 


•nee  of  a  self-consistent  real  frequency  shift  of  the  single  mode.  The  simultan¬ 
eous  constraints  imposed  by  F.qs  (34)  and  (35)  result  in  the  trapping  of  the  beam 
particles  and  the  corresponding  saturation  in  the  growth  of  the  single  mode. 

For(j*0,  it  is  evident  from  Eqs  (34)  and  (35)  that  a  steady  state  does  not 
exist.  Both,  the  total  momentum  and  the  total  energy  of  the  system  must  increase 
secularly  in  time  due  to  the  external  push.  There  are  two  conceptually  different 
states  capable  of  satisfying  Eqs  (34)  and  (35).  One  of  them  is  the  pure 
runaway  regime,  in  which  the  external  push  provided  by  the  DC  field  gives  rise 
to  an  unbounded  growth  in  beam  momentum,  i.e.,  dp^/dt^S.  The  other  state  is  the 
more  stA>tle  beam  clamping  regime,  in  which  the  beam  momentum  remains  essentially 
constant  while  the  wave  amplitude  increases  secularly  in  time,  i.e.,  dp^/dt'vS 
Of  course,  one  also  encounters  an  intermediate  state  in  which  only  a  small 
fraction  of  the  beam  remains  clamped,  thus  resulting  in  a  simultaneous  increase 

in  mechanical  and  wave  momentum. 

1 

The  attainment  of  a  given  state,  beam  or  wave  runaway,  is  determined  by  the 


strength  of  the  DC  electric  field.  If  F.o>>r.^.,  then  the  force  produced  by  the 


wave  can  r^Qt\cancel  the  external  push,  consequently  the  beam  runaway  regime  i: 


triggered.  However,  if  Eo<F..p,  then  the  push  of  the  external  field  can  be  reduced 


by  the  wave  electric  field.  As  is  discussed  in  Sec  II,  the  l>C  field  shifts  the 
position  of  the  trapped  particles  relative  to  the  bottom  of  the  well,  hence  the 
particles  continue  to  drive  the  wave  unstable.  Another  way  of  visualizing  this 
effect  is  to  say  that  the  DC  field  offsets  the  beam  recoil  associated  with  the 


saturation  process  for  F.o=0.  !;°r  I^Mi.p  it  is  expected  that  the  intermediate 


regime  should  set  in. 
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IV.  Clamping  Behavior 

To  investigate  the  dynamics  of  beam  clamping,  Eqs  (23)  and  (25)  are  solved 
numerically  for  a  collection  of  discretized  beam  particles,  typically  in  the 
range  N^300-500.  At  first  we  consider  the  undamped  (r^=0)  single  mode  problem 

in  Which  the  linearly  fastest  growing  mode  becomes  dominant  and  traps  the  beam 
electrons. 

In  Fig.  2,  the  smaller  curve,  labelled  Eo=0,  exhibits  the  time  evolution 
of  the  amplitude  of  the  single  mode  when  there  is  no  DC  electric  field  present. 

The  characteristic  long-lived  amplitude  oscillations  due  to  the  bouncing  motion  of 
the  trapped  particles  are  evident.  The  secularly  growing  solid  curve  in  Fig.  2 
is  the  result  obtained  by  turning  on  a  DC  field  of  strength  Eo=0.5E^,  at  the  time 
the  wave  reaches  the  first  peak,  i.e.,  at  the  saturation  time.  As  predicted  by 
the  spatially  averaged  formalism  and  required  by  the  conservation  laws,  the  wave 
exhibits  unbounded  growth  (wave  runaway)  while  the  momentum  of  the  beam  particles 
remains  clamped  (i.e.,  it  does  not  increase).  The  dashed  curve  in  Fig.  2  repre¬ 
sents  an  approximate  prediction  that  would  be  made  by  the  spatially  averaged 
formalism,  i.e., 

\.  »?  CO  +  6 IX-  Vm  1  (36, 

where,  TQ  refers  to  the  turn-on  time.  Of  course,  the  more  sophisticated  prediction 
of  the  spatially  averaged  theory  would  require  the  simultaneous  solution  of  3 
algebraic  equations,  as  is  done  in  Ref.  10.  Nevertheless,  the  simple  prediction 
of  Eq  (3o)  gives  a  good  account  of  the  average  secular  growth  induced  by  the  DC 
field.  However,  Eq  (36)  can  not  reproduce  the  oscillations  that  accompany  the 
secular  growth.  These  oscillations  are  associated  with  the  bouncing  motion  of 
the  trapped  particles  in  a  potential  well  which  grows  in  time,  as  described  by 
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the  analytic  solution  [Eqs  (20)  and  (22)]  of  the  model  discussed  in  Sec  II.  As 
expected.  Fig.  2  shows  that  the  frequency  of  the  amplitude  oscillations  increases, 
and  their  depth  disappears. 

Figure  3  displays  the  phase-space  patterns  at t ■35.0  for  the  two  cases  shown 
in  Fig.  2.  Note  that  the  velocity  scales  are  different  for  the  two  cases,  while 
the  spatial  scales  are  the  same.  In  Fig.  3  the  solid  curve  superimposed  on  the 
phase-space  pattern  is  the  instantaneous  spatial  dependence  of  the  force  produced 
by  the  wave.  The  major  difference  between  the  two  cases  is  that  for  EQ*0,  the 
trapped  particles  are  smeared  out  over  a  full  wavelength,  whereas  for  E^  O.SE^ 
the  particles  are  focused  and  exhibit  a  spiraling  motion  of  ever  decreasing 
spatial  extent.  It  should  be  stressed  that  the  average  momentum  of  these  particles 
does  not  increase  even  though  the  DC  field  is  acting  on  them.  In  addition,  the 
phase- space  for  Eq«  O.SE^,  shows  the  existence  of  a  secondary  population  of  run¬ 
aways  which  consists  of  less  than  10%  of  the  original  beam  particles. 

An  element  of  potential  weakness  in  the  spatially  averaged  theory  of  clamping 
is  the  need  to  specify  the  spatially  averaged  distribution  function  <f>  of  the 
beam,  as  it  interacts  with  the  growing  waves.  The  simplest  mathematical  choice 
is  used  in  that  formalism.  It  consists  simply  of  a  Lorentzian  whose  parameters 
change  self-consistently  with  the  conservation  laws.  The  reason  why  the  spatially 
averaged  formalism  is  so  successful  in  predicting  the  clamping  behavior  in  spite 
of  this  difficulty  is  that  the  actual  <f>  obtained  by  following  the  exact  particle 
orbits  has  the  essential  qualitative  features  of  the  Lorentzian  beam,  as  can  be 
seen  in  Fig.  4.  It  is  observed  in  this  figure  that  the  majority  of  the  particles 
form  a  clamped  beam  that  has  a  finite  thermal  spread  which  increases  slowly  in 
time.  Figure  4  also  shows  a  feature  which  is  completely  missed  by  the  spatially 
averaged  theory,  namely,  a  small  secondary  runaway  beam  which  can  not  be  stopped 
by  the  main  wave.  The  secondary  beam  can  interact  with  other  modes  having  faster 


-20- 


phase  velocities  and  nay  divide  itself  into  a  clamped  component  and  yet  another 

secondary  beam.  Therefore,  a  percolation  process  can  be  established  in  phase- 

space  that  produces  a  multiple  beam  distribution  function  of  the  type  observed 

22 

in  a  toroidal  strong  turbulence  experiment. 

The  enhanced  damping  effect  produced  by  reversing  the  direction  of  the  DC 
field  is  shown  in  Fig  Sa).  In  here,  Eoa-0.SEj.,  hence  the  beam  is  decelerated 
out  of  resonance  with  the  wave,  i.e.,  the  wave  switches  from  the  unstable  slow 
beam  root  to  the  damped  fast  component,  as  is  described  by  the  spatially  averaged 
theory.  The  originally  trapped  particles  envolve  into  a  decelerating'  beam  whose 
phase-space  structure  is  displayed  in  Fig  5b)  for  t*35.0. 

The  time  development  of  the  pure  beam  runaway  regime  is  exhibited  in  Fig. 

6a)  for  Eo=*4.0E^,;  the  Eq=0  is  included  for  comparison.  It  is  seen  that  the  DC 
field  initially  enhances  the  wave  amplitude  above  the  normal  saturation  level 
due  to  the  favorable  relative  displacement  of  the  bottom  of  the  potential  well 
described  previously.  However,  in  this  case  the  push  given  by  the  DC  field  is 
too  strong  and  cannot  be  compensated  by  the  enhancement  in  the  wave  amplitude. 

Accordingly,  the  trapped  particles  are  pushed  out  of  the  wave  troughs  and  become 
runaway  electrons  whose  phase-space  structure  is  shown  in  Fig  6b).  The  transition 
from  trapped  particles  to  runaways  encountered  in  this  problem  is  a  classical 
analog  of  the  process  of  ionizing  an  atom. 

The  clamping  behavior  exhibited  in  Fig.  2  has  been  obtained  by  turning  on 
the  DC  field  at  the  time  of  saturation  for  the  single  mode.  Therefore,  it  is  of 
interest  to  investigate  if  the  attainment  of  this  state  depends  crucially  on  this 
particular  choice  of  timing.  The  result  of  such  an  investigation  is  shown  in 
Fig.  7,  where  it  is  seen  that  the  clamping  state  characterized  by  the  secular 
increase  in  wave  amplitude  is  attained  over  a  wide  choice  of  initial  timing.  The 
timing  window  sampled  in  Fig.  7  ranges  from  the  pre-trapping  stage  (1)  through  » 


i 


-21- 


the  damping  (2-4)  und  rogrowth  (4-6)  cycles  of  the  dynamics.  The  principal 
effect  associated  with  the  choice  of  timing  is  the  natural  delay  of  the  onset 
of  the  secular  growth  stage,  which  is  manifested  in  Fig.  7  as  the  relative 
slippage  of  the  wave  amplitudes  for  the  different  cases  (1-6). 

The  sensitivity  of  the  clamping  state  to  the  amplitude  of  the  applied  DC 
electric  field  can  be  quantified  by  means  of  the  ratio  of  the  number  of  runaway 
electrons  created  NR  to  the  total  number  of  initial  beam  particles  N.  The  depen¬ 
dence  of  this  ratio  on  Eq/Ej  is  displayed  in  Fig.  8.  It  is  seen  that  for  E ^E^.cl.25 

the  runaways  constitute  less  than  10%  of  the  total  population  and  their  number  is 

not  sensitive  to  the  exact  value  of  E  .  The  reason  for  this  behavior  is  that  in 

o 

this  regime  a  significant  fraction  of  the  runaways  are  created  by  the  natural  spil¬ 
lage  of  the  trapped  particles  out  of  the  time  dependent  potential  wells;  a  process 
which  occurs  even  when  E  *0.  For  E  >1.25E_,  one  observes  a  dramatic  increase  in 
the  number  of  runaways,  eventually  attaining  the  100%  level  at  Eq=4.0Et.  This 
absolute  upper  limit  on  the  destruction  of  trapping  by  a  DC  field  is  roughly  a 
factor  of  2  larger  than  expected  for  a  wave  of  fixed  amplitude.  The  reason  for 
this  higher  threshold  is  that  the  wave  amplitude  grows  due  to  the  application  of 
the  DC  field,  hence  the  threshold  discussed  in  Sec  II  (i.e. ,  Eq=2E^.)  varies  as  a 
function  of  time  during  the  interaction.  The  net  integrated  result  is  essentially 
the  curve  shown  in  Fig  7. 

Since  the  clamped  state  is  insensitive  to  the  timing  of  the  DC  field  and  it 
exists  over  a  broad  range  of  values  Eq<4.0Et,  it  is  then  possible  to  modulate  the 
amplitude  of  the  wave  by  applying  DC  pulses  which  are  consistent  with  these  limits. 
An  example  of  such  amplitude  control  is  demonstrated  in  Fig  11,  where  one  applies 
two  square  pulses  of  magnitude  |Eo|-0.5Ej.,  but  having  opposite  polarity.  It  is 
seen  that  the  amplitude  can  be  increased,  held  constant,  and  decreased.  An  impor¬ 
tant  point  about  such  a  control,  Which  may  have  some  practical  application,  is 


that  the  response  of  the  wave  amplitude  is  essentially  linear  even  though  the 
dynamics  of  the  underlying  system  is  strongly  nonlinear. 

In  the  clamping  behavior  previously  discussed,  it  is  explicitly  assumed  that 
a  single  mode  dominates  the  wave  spectrum.  In  practice  this  situation  can  be  re¬ 
alized  if  the  system  evolves  from  an  extremely  low  noise  level,  or  if  the  spectrum 
is  initially  enhanced  by  modulating  the  beam  or  launching  a  test  wave  of  substan¬ 
tial  amplitude  and  with  a  wavenumber  corresponding  to  the  fastest  growing  linearly 
unstable  mode.  The  presence  of  many  modes  in  the  beam  instability  produces  a 
natural  smearing  of  the  phase-space  granulation  associated  with  particle  trapping. 
This  effect  can  be  viewed  dynamically  in  terms  of  a  non-stationary  potential  well 
(due  to  the  sum  over  all  modes)  which  gives  rise  to  the  stochastic  trapping  and 
detrapping  of  the  resonant  particles,  or  alternatively  as  arising  due  to  global 
phase  mixing  because  particles  trapped  in  different  troughs  bounce  with  different 
phases  and  bounce  frequencies.  In  either  interpretation  the  net  result  is  that 
the  cold  beam  develops  a  thermal  spread  rather  rapidly.  Therefore,  it  is  of  interest 
to  check  if  the  secular  wave  growth  triggered  by  the  application  of  the  DC  electric 
field  occurs  in  the  more  complicated  multi-wave  environment. 

The  role  played  by  modes  having  wavelengths  longer  than  that  of  the  fastest 
growing  mode  (i.e.,  k  <kQ)  is  shown  in  Fig  10.  This  figure  exhibits  the  time 
evolution  of  the  wave  amplitudes  of  a  system  containing  7  modes.  The  spectrum  con¬ 
sists  of  the  fastest  growing  mode  kQ,  and  6  lower  modes  having  a  mode  spacing 
Ak/kQ»0.07,  and  with  a  beam  density  corresponding  to  n/iOpB0.1.  In  Fig  10a)  the 

behavior  for  E  =0  is  shown.  It  is  observed  that  the  main  mode  is  still  dominant 
o 

throughout  the  entire  time  evolution.  However,  the  trapped  particle  amplitude 
oscillations  are  not  well  defined.  The  reason  for  this  behavior  is  seen  in  the 
phase-space  shown  in  Fig.  10  b),  where  a  high  degree  of  phase  mixing  is  observed. 
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The  effect  produced  by  the  application  of  a  DC  electric  field  (Eo>0.5Ej.)  to  this 
system  is  shorn  in  Fig.  11a).  It  is  found  that  the  characteristic  secular  growth 
in  wave  amplitude  indeed  occurs,  and  that  cleaT  trapping  oscillations  now  become 
evident.  Note  the  difference  in  amplitude  scales  between  Figs.  10a)  and  11a). 

As  expected  from  Fig.  Ha),  the  phase-space  associated  with  the  secular  growth. 

Fig.  lib),  shows  heavily  populated  clumps  and  a  negligible  numbeT  of  secondary 
runaways,  quite  similar  to  the  behavior  obtained  for  the  pure  single  mode  case. 

The  effect  produced  by  the  inclusion  of  waves  having  wavelengths  shorter  than 
that  of  the  main  mode  (i.e.,  k>kQ)  is  shown  in  Fig  12a).  The  spectrum  in  this 
system  contains  the  fastest  growing  mode  kQ,  and  two  upper  and  two  lower  sidebands 
with  mode  spacing  Ak/kQ=0.08,  and  fl/u^aO.l.  The  principal  feature  in  this  system 
is  that  the  main  mode  decays  right  after  reaching  saturation  and  triggers  the 
growth  of  an  upper  sideband  which  eventually  dominates  the  spectrum.  In  this 
case  the  characteristic  trapping  oscillations  are  absent,  as  is  expected  from  the 
corresponding  phase-space  shown  in  Fig  12b).  The  application  of  a  DC  electric 
field  (Eq=0. 7Ep)  to  this  system  also  produces  a  secular  growth  in  the  wave  amplitude, 
as  shown  in  Figure  13a).  Again,  note  the  change  in  scales  between  Figs  12a)  and 
13a).  The  degree  of  clump  formation  in  this  case  is  much  less  than  in  the  pure 
single  mode  problem  and  the  number  of  secondary  runaways  is  significantly  increased 
as  seen  in  the  phase-space  of  Fig  13b).  The  reason  for  this  behavior  is  that  in 
the  zero  order  system  (EQ*0),  particle  trapping  is  destroyed  by  the  sideband 
growth. 
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V.  Dynamic  BGK  Equilibrium 

It  has  been  shown  in  the  previous  section  that  in  the  absence  of  dissipation 
the  wave  amplitude  exhibits  unbounded  growth  for  DC  fields  such  that  Eo<4.0E». 
Within  the  restrictions  of  the  model  a  steady  state  cannot  be  attained.  Therefore, 
the  clamped  beam  particles  are  compressed  into  singular  charge  clumps.  However, 
when  a  finite  amount  of  wave  dissipation  is  present,  an  interesting  nonlinear 
steady  state  can  be  found  analytically.  It  consists  of  a  trapped  particle  equili¬ 
brium  of  the  type  associated  with  a  BGK  mode,  but  in  which  there  is  energy  and 
momentum  being  continuously  circulated  through  the  system. 

The  steady  state  flow  of  energy  and  momentum  can  be  deduced  from  Eqs  (27) 
and  (28)  for  the  single  mode  problem,  i.e., 

r- 1  =  6  "  *p  ,M'  (37) 


&_ 

AT 


9 


(38) 


where  r  is  the  damping  rate,  W  the  real  frequency  shift  and  A  the  wave  amplitude. 
The  steady  state  is  obtained  by  balancing  the  dissipation  of  energy  and  momentum 
due  to  the  wave  damping  against  the  source  of  momentum  and  energy  provided  by 
the  DC  field.  The  requirements  are 

a «(«/*■  >  w=p*  w 

These  conditions  aiise  from  global  properties  of  the  wave-particle  system  and  must 
be  consistent  with  the  motion  of  the  individual  particles  in  order  for  the  steady 
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state  to  be  realizable.  To  check  this  consistency  requirement,  write  the  equa¬ 
tion  of  motion  for  particle  j  in  real  form 

S  8  -  2  kb)  VliM  -  (40) 

dt*  4  • 

from  which  it  is  seen  that  dp^/dr-0  if  the  particle  position  is  given  by  (t)= 
e+Pbr.  In  here  6  is  a  constant  phase  factor  independent  of  the  particle  label 
j  (i.e.,  a  singular  charge  clump  is  formed),  and  must  satisfy  the  force  equili¬ 


brium  condition 


.5-  -  oae 
2  A 


which  has  a  physical  solution  provided fc/2A<l.  This  consistency  requirement, 
together  with  Eq  (39),  implies  that(€  r/2)^2<l  for  the  dynamic  BGK  equilibrium 
to  exist. 

It  must  then  be  verified  that  the  choice  of  a  singular  clump  indeed  gives 
rise  to  a  constant  wave  amplitude  and  a  constant  frequency  shift.  This  check 
can  be  obtained  by  analyzing  the  real  version  of  Eq  (25),  namely 


Ax  H 

AW. 

*  1 


The  choice  of  a  singular  clump  population  transforms  Eq  (42)  into 


-26- 


A.  ^  *  P  k  =  OCr*  0 

Jx 


which  using  the  dynamical  equilibrium  condition  of  Eq  (41),  together  with  the 
conservation  constraints  of  Eq  (39),  results  in 


A  k  *  PA  =.  ('2PA?V2A 
ax 


(45) 


thus  showing  that  dA/dt®0,  as  required. 

Similarly,  Eq  (43)  is  transformed  into 

Vyj  -  ©  /  A  (46) 

which  is  explicitly  time  independent  in  view  of  Eq  (45),  and  from  which  one  obtains 
the  self-consistent  frequency  shift  associated  with  the  dynamic  BGK  equilibrium 

wi=  -  (*r/£)l/z  £i-  (47) 

The  choice  of  the  negative  sign  in  the  square  root  of  Eq  (47)  is  determined  by 
the  physical  need  to  locate  the  clump  within  the  confining  trough  of  the  wave 
potential.  The  disregarded  positive  root  corresponds  to  a  potential  hill  instead 
of  a  valley. 

The  nonlinear  steady  state  analyzed  in  this  section  consists  of  a  singular 
charge  clump  which  is  located  inside  the  potential  trough  of  a  constant  amplitude 
wave  and  propagates  at  constant  speed  even  though  an  external  DC  electric  field 


-2 


! 

| 

acts  on  it.  This  is  an  example  of  a  nonlinear  dynamical  equilibrium  in  which 
the  trapped  beam  particles  play  the  role  of  an  intermediate  agent  which  transfers 
the  external  momentum  push  of  the  DC  electric  field  to  those  physical  mechanisms 
responsible  for  the  wave  damping.  The  physics  behind  the  continuous  excitation 
of  the  wave  is  that  the  constant  position  of  the  charge  clump  is  displaced  relative 
to  the  bottom  of  the  wave  trough  in  the  positive  energy  transfer  phase  of  the 
usual  trapped  particle  energy  transfer  cycle,  as  is  evident  from  Eq  (44). 

The  singular  charge  clump  associated  with  the  dynamic  BGK  equilibrium  can 
be  attained  only  in  the  limit  r-*0,  t-*»  (neglecting  space  charge  limits,  which 
eventually  become  important).  For  finite  values  of  r  and  t,  the  actual  clump 
has  a  finite  spatial  extent.  The  consequence  of  not  having  a  perfect  charge 
singularity  on  the  parameters  of  the  exact  dynamic  BGK  equilibrium  can  be 
assessed  by  expressing  the  orbit  of  the  j  trapped  particle  as 

*i=  \  c<») 

where  is  the  force  free  equilibrium  position  defined  by  cost  *  €/2A,  and 
Oj  is  the  finite  excursion  having  frequency  and  random  phase  0^. 

Considering  that  a  fraction  o<l  of  the  original  beam  particles  remain 
clamped  and  a  corresponding  amount  (1  -  a)  runs  away,  as  is  always  the  case, 
implies  that  the  self-consistent  frequency  and  amplitude  are  determined  essen- 


tially  by 

N  % 

(49) 

VI  K  a  £■  + 

N  1 

(50) 
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in  which  the  snail  (high  frequency)  modulation  due  to  the  oscillatory  orbits 
of  the  runaways,  given  in  Eq  (13),  is  neglected. 

Using  the  finite  series  expansion  in  terms  of  Bessel  functions  for  the 
expressions  in  Eqs  (49)  and  (50),  and  realizing  that  summing  over  j  is  equivalent 
to  averaging  over  the  random  phases  ^  leads  to 


PA  *  f 

(51) 

W  A  =  5  s**- 

(52) 

11 

•rM. 

T 

•A 

(53) 

where  JQ  is  the  Bessel  function  of  order  0  and  the  quantity  f  measures  the 
degree  of  focusing  of  the  clump,  e.g.,  as  a^-*0,  f*o. 

Solving  Eqs  (51)  and  (52)  leads  to  the  modified  parameters  foT  the  approx 
imate  BGK  equilibrium  due  to  a  finite  size  clump 

A  =  U  6 /*!’)"'  (54) 

IPe/iot'*  (S5) 

which  reduce  to  Eqs  (39)  and  (47)  in  the  limit  f»l.  From  Eq  (55)  one  finds 
the  more  restrictive  condition  for  the  existence  of  the  finite  size  BGK  equili 
brium  to  be  (r6/2) <£.  This  condition  should  be  used  in  conjunction  with  Fig  8 
to  delineate  the  parameter  space  ameanable  to  the  formation  of  a  dynamic  BGK 
equilibrium. 


To  check  on  the  existence  of  the  dynamic  BGK  equilibrium,  we  solve  Eqs  (23) 

and  (25)  as  in  Sec  IV,  but  now  proceed  to  include  the  effect  of  dating,  i.e., 

tyO.  Figure  14a)  shows  the  time  evolution  of  the  wave  amplitude  for  the  case 

r-0.02,  Eo«1.0  Ep.  It  is  found  now  that  the  secular  behavior,  exemplified  by  Fig. 

2,  is  stopped.  The  system  approaches  asymptotically  a  steady  state  whose  amplitude 

1/2 

is  slightly  smaller  than  the  prediction  of  Eq  (39),  i.e.,  (€/2r)  ■  4.2.  Sim¬ 

ultaneously  with  the  approach  to  this  steady  state,  one  observes  the  formation 
of  a  high  density  cluster  of  beam  particles,  a  dump,  whose  velocity  remains 
essentially  constant,  as  is  characteristic  of  the  dynamic  BGK  equilibrium,  and 
is  exhibited  in  the  phase-space  plot  of  Fig  14b). 

The  formation  of  a  highly  coherent  trapped  particle  state  in  the  presence  of 
finite  damping,  as  demonstrated  in  Fig.  14,  may  on  first  thought  appear  to  contra¬ 
dict  the  well  documented  fact  that  trapping  can  be  destroyed  due  to  the  effect  of 

0 

finite  damping,  as  studied  by  Dimonte  and  Malmberg.  As  shown  by  these  authors, 
the  destruction  of  trapping  is  not  due  to  the  obvious  dissipative  decrease  in  wave 
amplitude,  but  rather  it  is  caused  by  a  sudden  nonlinear  phase  shift  which  pulls 
the  particles  out  of  the  wave  troughs,  thus  causing  a  dramatic  damping  in  the  wave 
amplitude.  Figure  15  shows  the  existence  of  this  sudden  damping,  for  I*  =  0.075, 
occuring  just  after  the  second  trapped  particle  oscillation,  x<20.0.  The  phase- 
space  structure  corresponding  to  this  damping  is  shown  in  Fig.  16a),  where  it  is 
evident  that  strong  phase  mixing  has  occurred.  The  effect  produced  by  the  applica¬ 
tion  of  a  DC  electric  field  to  the  phase-space  structure  of  Fig  16a)  is  exhibited 
in  the  late  time  (t>80)  portion  of  Fig  15.  The  DC  field  causes  the  regrowth  of 
the  wave  and  asymptotically  it  approaches  a  dynamic  BGK  equilibrium.  The  regrowth 
of  the  wave  is  accomplished  by  the  reorganization  of  the  phase-space  structure  due 
to  the  push  provided  by  the  DC  field.  The  DC  push  positions  the  majority  of  the 
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trapped  particles  in  the  correct  phase  for  transferring  energy  to  the  wave,  as 
seen  in  the  phase-space  of  Fig  16b).  In  the  process  of  Moving  the  trapped  particles 
to  the  favorable  position  for  wave  growth  a  secondary  runaway  population  is  created. 

The  attainment  of  a  steady  state  due  to  finite  wave  dissipation  in  the  presence 
of  many  modes  is  also  possible,  as  shown  in  Fig.  17a)  for  the  parameters  Eo*0.7Ej., 
r-0.02.  The  system  contains  6  lower  modes  (kdcQ)  in  addition  to  the  main  mode,  and 
has  a  mode  spacing  Ak/k^O.07  with  fl/Wp«0.1,  as  In  the  case  considered  earlier  in 
Fig.  12.  Together  with  the  approach  to  an  asymptotic  constant  wave  amplitude,  one 
observes  in  the  corresponding  phase-space  picture  of  Fig.  17b)  that  spatially  focused 
charge  clumps  appear,  as  expected  from  the  dynamic  BGK  equilibrium  analysis. 

However,  due  to  the  presence  of  several  non-resonant  modes,  a  significant  fraction 
of  secondary  runaways  is  created.  Since  the  runaways  do  not  contribute  to 
driving  the  wave  amplitude,  it  is  found  that  the  asymptotic  wave  amplitude  attains 
a  level  A“3.0  which  is  lower  than  the  ideal  pure  single  mode  prediction  of  A=4.2. 

From  this  difference  it  can  be  estimated  through  Eq  (54)  that  the  focusing  parameter 
has  a  value  f“0. 56.  This  value,  although  lower  than  that  for  a  perfect  dynamic 
BGK  equilibrium  (f»l. 0) ,  is  found  to  satisfy  the  necessary  existence  criterion 
for  the  steady  state,  i.e.,  r(^2<f.  A  similar  result  is  also  obtained  when 
upper  sidebands  as  well  as  lower  sidebands  are  present.  However,  due  to  the 
destruction  of  particle  trapping  produced  by  the  upper  sidebands,  as  seen  in 
Fig.  12,  the  degree  of  clumping  is  significantly  reduced,  i.e.,  the  focusing 
parameter  attains  a  low  value  f*0.2. 


-31- 


VI.  Stability  of  Clump  Equilibrium 

Several  years  ago  considerable  theoretical  effort  was  devoted  to  the  study 

of  the  sideband  instability  associated  with  trapped  particle  BGK  equilibria.  The 

notivating  reason  for  the  intensive  theoretical  activity  was  the  experimental 

23 

observation  of  Wharton,  et  al.  that  sideband  signals  were  spontaneously  excited 

when  a  monochromatic  large  amplitude  electron  plasma  wave  was  externally  launched 

in  a  collionless  plasma.  One  of  the  more  interesting  models  proposed  to  explain 

24 

the  experimental  results  was  advanced  by  Kruer,  et  al.  It  consisted  of  the 

ad-hoc  assumption  that  the  trapped  particles  were  bunched  into  a  singular  clump 

located  exactly  at  the  bottom  of  the  potential  well  of  the  wave.  The  earliest 

analysis  of  this  model  by  Kruer,  et  al.  showed  that  sideband  waves  were  indeed 

unstable;  the  physics  behind  the  instability  being  essentially  associated  with 

the  linear  dynamics  of  a  spatially  periodic  cold  beam.  Although  a  computer 
25 

particle  simulation  of  the  problem  by  Kruer,  et  al.  showed  that  instead  of 
a  singular  clump,  a  large  hole  was  formed  at  the  bottom  of  the  well,  the  clump 
model  yielded  predictions  which  were  in  rough  agreement  with  the  observed  side¬ 
band  growth.  The  analysis  of  Kruer,  et  al.  was  subsequently  corrected  by  Goldman 
26 

and  Berk  to  include  the  important  effect  produced  by  the  requirement  of  self- 

consistency,  i.e.,  a  nonlinear  frequency  shift  must  always  accompany  the  formation 

of  the  charge  clumps.  The  nonlinear  frequency  used  by  Goldman  and  Berk  was  taken 

27 

from  an  early  calculation  by  Bohm  and  Gross  ,  and  although  it  provides  a  logical 
improvement  upon  the  model  of  Kruer,  et  al.,  it  is  not  the  correct  frequency 
shift5' ®  associated  with  a  large  amplitude  electron  plasma  wave.  These  difficulties 
simply  stem  from  the  fact  that  the  clump  model  is  not  the  correct  dynamical  descrip¬ 
tion  of  the  launched  wave  problem.  Nevertheless,  it  is  an  interesting  system 
to  investigate. 
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The  dynamic  BGK  equilibrium  discussed  in  Sec  V  has  the  singular  clump  nature 
of  the  model  postulated  by  Kruer,  et  al.,  but  does  not  suffer  from  the  lack  of 
self-consistency  and  dynamical  origin.  The  wave  amplitude,  clump  position,  and 
nonlinear  frequency  shift  are  all  analytically  known.  In  addition,  the  formation 
of  the  clump  is  dynamically  possible  due  to  the  spatial  compression  provided  by 
the  secular  growth  of  the  wave  potential.  Accordingly,  the  dynamic  BGK  equili¬ 
brium  is  the  natural  system  for  calculating  the  stability  of  clump  equilibria. 

The  stability  analysis  that  determines  the  growth  of  small  waves  of  amplitude 
wavenumber  k,  and  frequency  u,  can  be  performed  in  a  manner  analogous  to  the 
calculation  of  Kruer,  et  al.  The  result  is 

kv.)  -  k-*  (56) 

in  which  vq  refers  to  the  time  asymptotic  velocity  of  the  wave  (different  from 
vq),  and  u>b  is  the  bounce  frequency  in  the  presence  of  the  DC  electric  field 
given  by  Eq  (10). 

In  general,  the  stability  problem  involves  the  solution  of  an  infinite 
determinant  which  represents  the  coupling  of  all  the  Floquet  (Bloch)  states 
k+tkQ,  as  indicated  in  Eq  (56).  However,  in  the  limit  of  small  the  low¬ 

est  order  nontrivial  coupling  which  is  often  investigated  consists  of  retaining 
only  the  l»0,2  terms.  The  simplified  system  consists  of 


0 


(57) 
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where  e(l)«e(k,<i»),e(2)«e(k-2ko,w-2ky  Q),  and 

__  <ru» 

(u-  hi?.)1-  Okl 

The  expansion  of  Eq  (57)  leads  to  the  dispersion  relation 


which  has  the  characteristic  structure  of  a  parametric  instability  problem. 

Taking  into  consideration  that  u>“(Dp+6ii>,  |  6uj  <«»>p,  and  recalling  the  weak 
cold  beam  scaling  given  in  Eqs  (24)  and  (26)  reduces  Eq  (59)  to 


in  which  the  unknown  complex  frequency  has  been  scaled  to  and  the  k 

dependence  enters  through 


■T  «  (si) 


It  should  be  mentioned  that  in  arriving  at  Eq  (60),  the  background  plasma 
has  been  taken  to  be  cold,  hence  the  scaling  factor  0-Up (n^/2no) . 

The  factor  W  appearing  in  Eq  (61)  is  the  self-consistent  nonlinear  frequency 
shift  given  by 


VI  =  -  (»r«70''*  Cs/iOl'11 


(62) 
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and  the  scaled  bounce  frequency  is 


Wb  =  (2  <rt/p)/4  [ \ -  (P6/2<r)]/4[v-  (63) 

The  principal  features  of  the  dispersion  relation,  Eq  (60),  can  be  obtained 
by  considering  the  limit  |w|-K),  \>l.  Physically,  this  regime  corresponds  to  weak 
damping  and  large  potential  wells,  as  is  required  for  the  formation  of  a  dynamic 
BGK  equilibrium.  Neglecting  the  explicit  dependence  on  r  and  W  in  the  left  side 
of  Eq  (60)  yields 


q-iw\ _ _ 

which  shows  that  as  |w|-*0  the  regions  of  instability  are  centered  around  r=sNb. 
Considering  each  of  these  regions  to  be  well  separated  from  each  other  results 
in 


*  1  <r,w' 


(65) 


where  r*s±±N^.  one  can  now  define  the  useful  scaled  variable 

(*fc*S*y  (66) 

to  express  Eq  (65)  in  terms  of  a  dispersion  relation  which  is  identical  to  that 
of  the  weak  cold  beam-plasma  instability 
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t\  (67, 

with  tc±-s±/(|w|o/2K|})1^3.  As  is  well  known,  the  maximum  growth  rate  in  Eq  (67) 
occurs  for  ic±«0,  and  has  growing  roots  given  by  exp  (i2ir/3)  for  and 

and  z_«exp(iir/3)  for  r^-W^.  This  implies  that  the  maximum  growth  rate  of  the 
instability  is 

and  corresponds  to  two  growing  sidebands 

=  v.  (  * * 

*  \  AVl/U)^  / 

(69) 


centered  around  the  main  carrier  mode  kQ,  and  each  having  a  spectral  width 


(  It  vO 


(70) 


Using  the  self-consistent  values  for  W  and  given  in  Eqs  (62)  and  (63), 
yields  the  parametric  dependence  of  the  growth  rate  y  from  Eq  (68),  namely 

t/Uf» 


(71) 
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This  result  shows  that  the  sideband  instability  driven  by  the  charge  clumps 
associated  with  the  dynamic  Bilk  equilibrium  disappears  asymptotically  as  r/$tO. 

For  finite  values  of  T/  ,  two  growing  sidebands  appear  at  the  bounce  wavenumber 

The  existence  of  these  two  unstable  sidebands  is  illustrated  in  Fig  18,  which 
is  obtained  by  the  exact  numerical  solution  of  the  dispersion  relation,  Eq  (60), 
for  the  parameters  I’=0.02,  $  =1 . 0,o=l . The  location  of  the  sidebands  coincides 
with  the  prediction  (k^)+= (1 . 3)ko>  (k^)  =(0.7)ko.  However,  the  maximum  growth 
rate  obtained  (c0.2S)  is  slightly  smaller  than  that  predicted  by  Eq  (68)  (i.e., 
,'/0.31).  The  difference  can  be  traced  to  the  fact  that  the  term  2|w|=0.4,  hence 
it  is  not  begligible  in  the  left  side  of  Eq  (60),  as  is  assumed  in  generating 
Eq  {**#),  which  is  strictly  correct  in  the  limit  | W | ->0 . 
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VII.  Conclusions 

This  study  has  demonstrated  that  the  ’phenomenon  of  beam  damping,  predicted 
by  the  spatially  averaged  formalism,  indeed  occurs  when  the  exact  microscopic 
dynamics  of  the  beam  particles  are  calculated  self -consistently.  The  clamping 
in  beam  energy  produces  an  unbounded  increase  in  wave  amplitude.  Therefore,  the 
spatial  structure  associated  with  beam  clamping  (and  averaged  out  in  the  previous 
formalism)  corresponds  to  trapped  particles  which  are  strongly  focused  in  space 
due  to  the  large  growing  potential  wells.  This  process  can  have  practical  applica 
tions  in  the  areas  of  efficiency  enhancement  of  travelling  wave  tube  devices  and 
future  free  electron  lasers. 

The  limits  on  the  amplitude  of  the  DC  fields  which  can  be  applied  before  run¬ 
away  behavior  sets  in  have  been  found.  A  useful  insensitivity  of  the  wave  ampli¬ 
tude  enhancement  on  the  time  of  application  of  the  DC  field  has  been  demonstrated. 

In  the  presence  of  wave  dissipation  a  new  nonlinear  state  has  been  found 
analytically,  and  is  observed  to  be  approached  in  the  time  asymptotic  limit  by 
the  numerical  solutions  of  the  problem.  This  state  is  a  dynamic  BKG  equilibrium 
consisting  of  a  singular  charge  clump  which  propagates  through  the  medium  at 
constant  speed  even  though  a  DC  electric  field  acts  on  it.  The  clump  plays  the 
role  of  an  intermediate  agent  which  transfers  the  external  push  to  the  sources 
responsible  for  the  wave  dissipation. 

The  dynamic  BGK  equilibrium  found  in  this  work  provides  the  ideal  analytical 
model  to  investigate  the  sideband  instability  triggered  by  trapped  clumps,  i.e., 
all  its  parameters  are  known  analytically  and  it  can  be  generated  dynamically. 

The  growth  rate  for  this  instability  has  been  found. 

In  summary,  a  variety  of  interesting  effects  have  been  shown  to  arise  when 
a  DC  electric  field  is  applied  to  a  system  containing  trapped  particles.  Many 
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of  these  effects  appear  to  be  readily  aaenable  to  experimental  observation,  and 
some  may  find  a  practical  application. 
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Fig.  1: 


Fig.  2 


Fig.  3 


Fig.  4 


Fig.  5 


Figure  Captions 

Spatial  Dependence  of  the  effective  trapping  well  in  the  presence  of  a 
DC  field  for  the  cases  Eq/A-0,  0.2,  0.5,  0.9. 

Secular  growth  of  wave  amplitude  produced  by  a  DC  field  Eo*0.5E_.  The 
Eo«0  curve  is  the  saturated  beam  instability  in  the  absence  of  the  DC 
field.  The  dashed  curve  is  the  approximate  prediction  based  on  con¬ 
servation  of  momentum  alone. 

Phase-space  at  t*35.0  for  a)  E  *0.5E_  and  b)  E  *0.  The  solid  curves 

o  T  o 

represent  the  instantaneous  spatial  dependence  of  the  wave  force. 

Spatially  averaged  velocity  distribution  function  in  the  clamped  stage 
(t»40.0)  for  EQ-0.5  E^.  Note  the  appearance  of  a  small  secondary 
runaway  component. 

Enhanced  wave  damping  produced  by  a  decelerating  electric  field, 
a)  Time  evolution  of  wave  amplitude  for  Eq»0  and  Eq»-0. 5E^..  b) 
Phase-space  at  t-35.0  for  Eq»-0. SE^. 


Fig.  6:  Runaway  beam  regime,  a)  Time  evolution  of  wave  amplitude  for 
Eq«0  and  E^.O&p  b)  Phase-space  att'35.0. 
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Fig.  7:  Dependence  of  secular  growth  of  wave  amplitude  on  turn-on  time 

of  DC  field.  Inset  shows  turn-on  time  relative  to  E  «0  evolution. 

0 

Labels  on  secularly  growing  curves  correspond  to  inset  timing  marks. 

Fig.  8:  Dependence  of  secondary  runaway  production  on  strength  of  DC  field. 

Fig.  9:  External  control  of  wave  amplitude  by  two  square- wave  pulses  of 
opposite  polarity,  E  ■±0.5ET. 

Fig.  10:  Effect  of  lower  sidebands  (k<kQ)  for  Eo*0.  a)  Time  evolution  of 
wave  amplitude  for  7  modes,  Ak/kQ«0.07.  b)  Phase-space  at 
x<*28.0. 

Fig.  11:  Induced  secular  growth  in  the  presence  of  lower  sidebands  (k<kQ) 
for  Eq*0. 5E^>.  a)  Time  evolution  of  wave  amplitude  for  7  modes, 
Ak/kQ*0.07.  b)  Phase-space  at  x  **26.0.  To  be  compared  with  Fig.  10. 

Fig.  12:  Effect  of  lower  (k<kQ)  and  upper  (k>kQ)  sidebands  for  EQ«0. 

a)  Time  evolution  of  wave  amplitudes  for  5  modes  (two  upper 
and  two  lower  sidebands),  Ak/kQ»0.08.  b)  Phase-space  at  t  <*26.0. 

Fig.  13:  Induced  secular  growth  in  the  presence  of  lower  (k<ko)  and  upper 
(k>ko)  sidebands  for  Eq-0.7Et.  a)  Time  evolution  of  wave 
amplitudes  for  5  modes  (two  upper  and  two  lower  sidebands), 
Ak/kQ**0.08.  b)  Phase-space  at  x"0.02.  To  be  compared  with 
Fig.  12. 
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Fig.  14:  Tine  asymptotic  approach  to  a  dynamic  BGK  equilibrium  in  the  presence 

of  dissipation,  E  «0.7E].,  r-0.02.  a)  Time  evolution  of  wave  amplitude. 
Theoretical  prediction  of  steady  state  amplitude  is  A*4.2.  b)  Phase- 
space  at  t*26.0. 

Fig.  IS:  DC  field  reversal  of  the  destruction  of  trapping  oscillations  caused 
by  finite  dissipation.  r*0. 075,  Eq*0.  SE^,.  Sudden  destruction  occurs 
at  t“18.0,  Eq  is  applied  at  t"28.0. 

Fig.  16:  Phase-space  structures  for  a  finite  damping  case  r=0.075  a)  Subsequent 
to  sudden  damping,  t«22.0,  Eq=0.  b)  After  DC  electric  field  is  applied, 
Eo»0.5E^.,  t*49.0. 

Fig.  17:  Time  asymptotic  approach  to  dynamic  BGK  equilibrium  in  the  presence 
of  lower  sidebands  (k<kQ),  Eo«0.07E^,  r»0.02.  a)  Time  evolution 
of  wave  amplitudes  for  7  modes,  Ak/ko»0.07.  b)  Phase-space  at 
t.36.0. 

Fig.  18:  Scaled  growth  rate  Yj  and  frequency  shift  YR  for  unstable  sidebands 

triggered  by  dynamic  BGK  equilibrium  defined  by  E  ■l.OBj.,  r»0.02,  a=1.0. 
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